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<N : Abstract 

<D , We examine the B — > 6K decays within the framework of SUGRA models making use of the 

00 . 

qq | improved QCD factorization method of Beneke et al. which allows calculations of non-factorizable 

(T) ■ contributions. All other experimental constraints {B — > X^, neutron and electron electric dipole 

> : 

CN . moments, dark matter constraints, etc.) are imposed. We calculate the CP violating parameters 

m ■ 

S(j>K S i C<t>K s an d -A^k^ as wen as the branching ratios (BR) of B° and B^, Br[B — > <ftK]. We 
find for the Standard Model(SM) and mSUGRA it is not possible to account for the observed 
2.7a deviation between S^Kg and Sj/^x s . In general the BRs are also in 3a disagreement with 
experiment, except in the parameter region where the weak annihilation terms dominate the decay 
(and hence where the theory is least reliable). Thus if future data confirm the current numbers, 



this would represent the first significant breakdown of both the SM and mSUGRA. We show then 



that adding a SUGRA non-universal A soft breaking left-right term mixing the second and third 
generations in either the down or up quark sector, all data can be accommodated for a wide range 
of parameters. The full 6x6 quark mass matrices are used and the SUSY contributions calculated 
without approximation. 
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I. INTRODUCTION 



Rare decay modes of the B meson are important places to test the Standard Model (SM) 
and to look for new physics. However, large theoretical uncertainties in the calculations 
of exclusive non-leptonic B decays make it difficult to extract useful information from ex- 
perimental data. Nevertheless, CP asymmetries of neutral B meson decays into final CP 
eigenstates, i.e., B — ► <pKs and B — > J/^Ks, are uniquely clean in their theoretical inter- 
pretations. Among these decay modes, B° — > <pKs is induced only at the one loop level in 
the SM and hence is a very promising mode to see the effects of new physics. In the SM, it 
is predicted that the CP asymmetries of B° — > <fiKs and B — > J/^Ks should measure the 

^ible 0(X 2 ) difference On the other hand, the BaBar and Belle 

11. 



with negligible 0(X 2 ) difference [lj|. On t 

2I Ely, 111 show a 2.7a disagreement between S^Ks anc ^ ^J/9K S U 



measurements 



Sj/v Ks = 0.734 ±0.055, 

= -0.38 ±0.41. (1) 



while Sj/^Ks = sin2 / 9j/^,^- s (which is a tree level process) is in excellent agreement with 
Buras' SM evaluation, sin 2(3 = 0.715lo;oi from the CKM matrix Q]. In addition, the 
branching ratios (BRs) and the direct CP asymmetries of both the charged and neutral 
modes of B — > <pK have also been measured 



asymmetric 

ffflflfr 

Br[B° -> <f>K s ] = (8.0 ± 1.3) x 10~ 6 , 



Br[B + -> 4>K + ] = (10.9 ± 1.0) x 10~ 6 , (2) 
C+ Ka = -0.19 ±0.30, 
A CP (B + ^(f)K + ) = (3.9 ±8.8 ±1.1)%. (3) 

In general, any model should explain all these data. In particular, the relatively small 
uncertainties in the BRs of B + — > <pK + and B° — > <pKs need to be considered in the 
analysis since they are highly correlated and both are based on the b — > s transition. In the 



1 After submitting this work new data from Belle p gave a value of S^Ks = —0.96 ± 0.5^q'5i (a 3.5<r 
deviation from the Standard Model) and preliminary analysis of new data from BaBar gave S^k s = 
+0.45 ± 0.43 ± 0.07. Belle and BaBar would then disagree by 2.\u and if one averages the new values one 



obtains [JJ S^Ks — —0.15 ± 0.33 which is again 2.7a from the Standard Model- 
2 Our average of Br[_B + — > <pK + ] only includes BaBar and Belle since CLEO [9j is 2.3ct away. Br[_B + 
<j)K+] would become (9.4 ± 0.9) x 1(T 6 if CLEO is included. 
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SM, Acp{B + — > cj)K + ) is small and agrees with (jSJ). So this direct CP asymmetry result 
plays an important role in constraining the new physics contribution which might explain 
the discrepancy between S^k s and Sj/q,K s - 

The dis crep anc y b etween S^ wc a nd Su$k s nas been discussed in some recent works jicl 



l22j in the framework of SUSY models, especially in 
the minimal supersymmetric standard model (MSSM) with the mass insertion method [23]. 
Although these works can provide useful constraints on certain off-diagonal terms of squark 
mass matrices at low energy, we find that it is very interesting to investigate this problem 
in the context of grand unified theory (GUT) models. Since GUT models explain a number 
of phenomena at low energy by a few well motivated parameters at the GUT scale, various 
experimental measurements get correlated in this framework. Among supersymmetric GUT 
models, the R-parity conserved SUGRA model is one of the most favored models since it 
provides a natural explanation of the dark matter problem. The minimal SUGRA model 
(mSUGRA) with R-parity conservation has been investigated extensively because of its 
predictive power that comes from the fact that it depends on only a few new parameters. 
Unlike the MSSM, which is hard to be constrained due to its more than 100 new parameters 
(including 43 CP violating phases), the parameter space of the minimum SUGRA model 
has 4 parameters and 4 phases. This parameter space has several experimental constraints, 
i.e. b — > s + 7, neutron and electron electric dipole moments (EDM), LEP bounds and 

nn 

relic density measurements |24l. 1251] etc.. In this paper, we examine the observed discrepancy 
between S$k s and Sj/^k s in the context of SUGRA models including mSUGRA and models 
with non-universalities. We consider all relevant experimental constraints in our calculation. 

The calculation of B — > <f>K decays involve the evaluation of the matrix elements of related 
operators in the effective Hamiltonian, which is the most difficult part in this calculation. 
However, the newly developed QCD improved factorization (BBNS approach) provides 
a systematic way to calculate the matrix elements of a large class of B decays with significant 
improvements over the old factorization approach (naive factorization). It allows a QCD 
calculation of "non-factorizable" contributions and model independent predictions for strong 
phases which are important in the theoretical evaluation of the direct CP asymmetries of B 
decays, e.g. for B~ — > (f)K~ , whose experimental result is given in (J3J). We adopt the QCD 



improved factorization in our B — > <pK calculations. Recently Du et al. 



27. 



28|, |2S! have 



published an improved calculation of B — > PV decays. We followed here their calculational 
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(a) (b) 

FIG. 1: Hard spectator scattering diagram (a) and weak annihilation diagram (b). In (a) the gluon 
can connect the spectator with either <j) quark and in (b) the gluon can originate from any B quark 
or K quark. 



techniques |22[ which are based on the original work [26] of Beneke, et al. While the 
BBNS approach is an important advance in calculating B decays, it is not completely model 
independent. In the BBNS approach the hard gluon (H) and annihilation (A) diagrams (see 
Fig^) contain infrared divergences which are parameterized by an amplitude Ph,a (with 
Ph,a < 1) an d a phase 4>h,a- (More details can be found in |26j and 27].) If the effects 
of these terms are small, the theoretical predictions are well defined. However, if these 
terms are large or dominant, the theory becomes suspect. We will see below that S^k s is 
essentially independent of the infrared divergent terms, though the branching ratios can 
become sensitive to pa and 4>a- 

In this work, we first examine the mSUGRA model which is universal at the GUT scale 
and then consider non-universal terms. Although non-universalities may cause serious prob- 
lems in some flavor changing processes, e.g. K° — K° mixing and b — > s + 7, we will show 
that the off-diagonal terms in the A parameter soft-breaking terms can satisfy all experi- 
mental data including the 2.7a deviation between S^Ks an d Sj/^k s - We calculate the BR 
of B~ — > <pK~ which is highly correlated with B° — > <j)K° and calculate the CP asymmetry 
of this mode in the allowed parameter space. We also calculate the CP asymmetry of the 
b — > S7 decay mode. 

This paper is organized as following: In Sec. 2, we give a brief description of the CP 
asymmetry of B° — * (j)Ks and the QCD factorization technique used in this paper. Then 
we discuss the SUGRA model and its contributions to B° — > 4>Ks in Sec. 3. Some detailed 
discussion on experimental constraints implemented in our analysis are given in Sec. 4 and 
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the Standard Model predictions are discussed in Sec. 5. Sec. 6 is devoted to our results for 
the mSUGRA model, after which we proceed to the models with non-universalities in Sec. 7 
and we give conclusions in the last section. 



II. CP ASYMMETRY OF B <pK DECAYS 

The time dependent CP asymmetry of B — > <pKs is described by: 

TCBphysW - <t>K S ) - T(B phjs (t) -> <j>K & 



r(fiphysW - <PKs) + r(B° phys (t) -> <j>K s ) 

-C^Kg cos(Am B t) + S^Ks sm(Am B t) (4) 



where S^Ks and C$k s are given by 

2 ImA^ s 1 — |A^ S ' 2 

; |2 ' "W*s - TTT\ 

and A<^ s can be written in terms of decay amplitudes 



- TZTT\ 12 ' - T— jT [2 > W 

t + Aak-J 1 + A</- 



" ~ 6 .4(5° _ 0/^) (6) 

In our analysis, we find that the SUSY contributions to B,i — Bj mixing are small, and so 
from now on we will use the standard definition for (3: 

Within the SM, sin 2(3 can be measured by Sj/^k s - The current experimental result is given 
in Eq. (0). Since B° — > J/^Ks decay is dominated by the SM tree level contribution, we 
expect that in our analysis the new physics will not affect the SM prediction for sin 2/3 from 
B — > J/^Ks- As a consequence, we further assume that the current SM fit for the CKM 
matrix will not be affected by models discussed in this paper. 
The CP asymmetry of charged B — > (pK decay is defined as 

A r(£T -> $K-) - T(B+ -> <j>K+) = |A^ T | 2 ~1 

***** " r(B- -> + r(5+ -> 0K+) |A^ T | 2 + i u 



where 
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From the above discussion, it is clear that our theoretical predictions for the experimental 
observables, e.g. S^Ks, C<i>k s and A^kt, depend on the evaluation of decay amplitudes where 
the effective Hamiltonian plays an important role. The Effective Hamiltonian for B — > <pK 
in the SM is Q: 



V ^ q=u, c 



10 

c 1 o q 1 + c 2 o q 2 + Y, c MOM + c 7l o 7l + c 8g o 8g +h. c . (10) 

k=3 

where the Wilson coefficients Ci(fi) can be obtained by running the RGE from the weak 
scale down to scale fi. The definitions of the operator Oj's in the SM can be found in 
2(ij | . The SUSY contributions will bring in new operators Oj's which can be obtained by 
changing L «-> R in the SM operators. We use Ci to denote the Wilson coefficient of Oj. 
The decay amplitude of B — > M\M 2 can be expressed in terms of the matrix elements of 
Oj's, {M\M 2 \Oi\B) . We evaluate these matrix elements in the QCD improved factorization 
technique. The necessary expressions can be found in |2fil l^. 
Using the above Hamiltonian the amplitude of B — > <pK is: 

A{B -> <j)K) = A f (B -> <j)K) + A a (B -> <f>K) (11) 



where A^ is factorized amplitudes which can be written as |2 

G p 



V ^ p=u, c i 



(12) 



and v4. a is the weak annihilation decay amplitudes |27l |: 



A a (B - 0JO = E W' ( 13 ) 

The matrix elements (0iC|Oj|S)/ in Eq.()12|) are the factorized hadronic matrix elements 
[3^. Oj's and 6j's contain the Wilson coefficients. Explicit expressions for them, as well as 
for A a [B -> 0J-ST), can be found in Q and (27 1. 

In our discussion, the dominant SUSY contributions occur through and 0$ g and we 
calculate these new SUSY contributions in the SUGRA framework. 



III. SUGRA MODELS 

The SUGRA model at the GUT scale can be described by its superpotential and soft- 
breaking terms: 

W = Y U QH 2 U + Y D QH 1 D + Y L LH l E + ^H X E 2 
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i \ Z a 

(A U QH 2 U + A D QH 1 D + A L LH 1 E^j + h.c.) . (14) 

Here Q, L are the left handed quark and lepton doublets, U, D and E are the right handed 
up, down and lepton singlets and Hi t2 are the Higgs doublets. In the minimal picture, the 
mSUGRA model contains a universal scalar mass m , a universal gaugino mass mm and 
the universal cubic scalar A terms: 

2 --i m a =m 1/2j A U ' D,L = A Y U ' D ' L . (15) 



m„- = m 



This model contains four free parameters and a sign: m , mi/ 2 , A , tan/? = (H 2 )/(Hi) and 
the sign of \i. 

However, the parameters m^, fi and A can be complex and their phases can be 0(1). 
In order to accommodate the experimental bounds on the electron and neutron electric 
dipole moments (EDMs) without fine tuning phases we extend mSUGRA by allowing the 
gaugino masses at Ma to have arbitrary phases. This model has been extensively studied 



in the literature 

rrii 



EL 



32 



331 ] . Thus the SUSY parameters with phases at the GUT scale are 
mi/^e 1 ^ i=l,2,3 (the gaugino masses are for the U(l), SU{2) and SU(3) groups), 
Aq = \Ao\e lotA and /i = e*^^ . However, we can set one of the gaugino phases to zero 
and we choose <p 2 = 0. Therefore, we are left with four phases. The EDMs of the electron 



31 



32, 331. In our 



and neutron can now allow the existence of large phases in the theory 
calculation, we use 0(1) phases but calculate the EDMs to make sure that current bounds 
are satisfied. 

We evolve the above parameters from the GUT scale down to the weak scale using full 
matrix RGEs. Since the b — > s transition is a generation mixing process, it is necessary to 
use the full 6x6 matrix form of squark mass matrices in the calculation. We perform the 
calculation of SUSY contributions without any approximation. 

We also include the one loop correction to bottom quark mass from SUSY Q, 

which is 

important in the calculation of SUSY contributions to the Wilson coefficients of the operator 
07 7 and Os g and consequently affects the calculations of B — >• X s ^ and B — > <f>K decays. 
We now discuss the experimental constraints in the next section. 
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IV. PARAMETER SPACE AND EXPERIMENTAL BOUNDS 



In this section we review all the experimental constraints considered in our analysis and 
briefly discuss their effects and importance. 



A. B - X sl 

We use a relatively broad range for the branching ratio of B — > X s , y to take into 
account the uncertainty in the theoretical calculation of B — > X s ^ (±0.3 x 10 -4 ): 

2.2 x 10~ 4 < Br(fi -> X sl ) < 4.5 x 10" 4 . (16) 



The SM prediction for the Br[S — > X s ^\ is very close to the measured value |3j|, so the 
b — > s transition in any new physics is strongly constrained. Since the B — > (j)K decay also 
depends on the b — > s transition, the Br [B — > X^] constraint needs to be implemented in 
any analysis of the B — > <pK decay. Besides the BR of B — > X s 7, we also consider in this 
work the direct CP asymmetry of B — > X s ^ for the existing of CP violating phases. The 
experimental measurement from CLEO gives 37| : 

A^s+ 7 = (-0.079 ± 0.108 ± 0.022)(1.0 ± 0.030) (17) 

or at 90% confidence level, -0.27 < A b ^ s+1 < +0.10. 



B. The relic density 



The recent WMAP result gives 




n CDM h 2 = 0.1126t°;{™. (18) 

We implement this bound at the 2a level in our calculation:0.094 < Vth 2 < 0.129. We also 
notice that when non- universal terms are present, new annihilation channels may arise and 
they are different from the usual mSUGRA f — x° co-annihilation channel. 
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C. K° -K° mixing 

It has been shown that K° — K° mixing can significantly constrain certain flavor changing 
sources in SUSY models 3|. The current experimental bound for AMk is 4(3 ]: 

AM K = 3.490 ± 0.006 x 10~ 12 MeV (19) 

In non-universal SUGRA models, even if non-universal terms between the first two genera- 
tions are not present at the GUT scale, AMk can become large. This happens because the 
degeneracy between the first two generations can get broken by other non-universal terms 
via the RGEs. For example, the A 32 terms in the trilinear coupling matrices A in Eq. (jl4|) 
can give rise to new contribution to the m\ 2 term via 8m\ 2 oc A2 3 A^2Log[MGUT I M wea k\. 
Therefore, it is important to pay attention to AMk even when there is no apparent direct 
source producing large contributions to AMk- 

D. Neutron and electron electric dipole moments 

Neutron and electron electric dipole moments (EDM) can arise in any model with new 
CP violating phases. In SUSY models, an electron EDM arises from diagrams involving 
intermediate chargino-sneutrino states and intermediate neutralino-selectron states (for more 
details, see The current experimental bounds on neutron and electron EDMs 



are 

d n < 6.3 x 10~ 26 e cm, d e < 0.21 x 10~ 26 e cm (20) 

There are other important phenomenological constraints considered, e.g. bounds on 
masses of SUSY particles and the lightest Higgs (m/j > 114 GeV). 

V. B — >■ <fiK DECAYS IN THE STANDARD MODEL 

We first discuss B — > <pK decays in the SM. The largest theoretical uncertainties in this 
calculation come from weak annihilation diagrams which mostly depend on the divergent 
end-point integrals Xa parameterized in the form j3, 1^2 1 

X A = (1 + p A e i<t>A ) In A h = Aqcd, Pa < 1. (21) 
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Hard spectator processes contain similar integrals Xh which are parameterized in the same 
way. However, uncertainties from the hard spectator calculation are much smaller than 
those from the weak annihilation for this decay, so we will mainly concentrate on the later. 
These weak annihilation contributions depend also on the strange quark mass, m s , through 
the chirally enhanced factor n x : 



K 



2m 2 K 



m b (m s + m q ) 



(22) 



where m q is or m u . 

In Fig|21we show the dependence of the branching ratio of B~ — > <pK~ mode on (pA and 
m s for pA = 1. Fig|3] shows the dependence of the BR on pa- Similar graphs can be obtained 
for B° — > <j)K°. Since m the SM (Eq.(jHJ)) is small, we can compare Br[5 _ — > (j)K~ 

with the experimental measurement of Br[i? + — > <pK + ] given in Eq.(J2J). Before we discuss 
the graphs, we first list our parameters: pn = 1, 4>h — —68°, fg = 180 MeV, = 233 MeV, 
f K = 160 MeV and F BK = 0.34. The CKM matrix elements can be obtained through the 



Wolfenstein parameterization |4l| with A = 0.819, A = 0.2237, p = 0.224 and r] = 0.324. The 
integral Jo ($b (£)/£) ^£ = m B/^B where is the B meson light-cone distribution amplitude 
is parameterized by A# = (0.35 ± 0.15) GeV. For ji = 2.5 GeV we use A# = 0.2 GeV, and 
for p = nib we use Xb = 0.47 GeV. In addition, we always use asymptotic forms of the 



meson light-cone distribution amplitudes j26j,|27(. If not mentioned, we will use the above 
parameters in later calculations. 

In both figures, we give results for two different scales and two different m s values, i.e., 
p = nib by solid lines (m s (2GeV) = 96 MeV) and the dot-dashed lines (m s (2GeV) = 
150 MeV) and p. = 2.5 GeV by dashed lines (m s (2GeV) = 96 MeV). One can see that the 
scale dependence is not significant. The straight lines correspond to the branching ratios 
neglecting the weak annihilation contribution. Comparing Figs 121 and El we see that a large 
branching ratio comparable to the experimental value is obtained in the region pa — 1 
and 4>a — 0(or 2ir). However, in this region the weak annihilation diagrams dominate the 
branching ratio and thus the theory is most suspect. In the remaining part of the parameter 
space, where the weak annihilation effects are small and the theory is presumably reliable, 
the SM prediction of the branching ratio is small, about 3a below the experimental value. 
We conclude therefore that where the theory is reliable the SM is in significant disagreement 
with the experimental value of the Br[5 + — > 4>K + ], and in order to obtain a SM value in 
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14 l— r 

12 - 




1 2 3 4 5 6 

FIG. 2: Branching ratio of B~ — > at = 1. The solid curve corresponds to fi = nib, dashed 
curve for fx = 2.5 GeV with m s (2GeV) = 96MeV and the dot-dashed curve for \i = mt, with 
m s (2GeV) = 150 MeV. The two straight lines correspond to the cases without weak annihilation. 




4 ■ ~" - 
2 ■ 

0.2 0.4 0.6 0.8 1 

FIG. 3: Branching ratio of B~ — > <pK~ at <f>A = 0. The solid curve corresponds to /i = m^, dashed 
curve for \i = 2.5 GeV with m s (2GeV) = 96 MeV and the dot-dashed curve for \i = m& with 
m s (2GeV) = 150 MeV. The two straight lines correspond to the cases without weak annihilation. 
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accord with the experiment one must use parameters where the theory is least reliable. A 
similar result holds for the Br [B° — > (f)K s ]. Here theory predicts a branching ratio about 
10% smaller than for B + — > <j)K + (in accord with the experimental values of Eq.(j2J)) but the 
SM can achieve this only in the region where the weak annihilation processes dominate. 

The dot-dashed line, in FigJ21 corresponds to a larger value of m s and we see that the 
BR is very sensitive to m s only in the large annihilation region. The region with sufficiently 
large annihilation to accommodate the data decreases as m s increases, since the annihilation 
amplitude then decreases, as can be seen from Eq. (j2"2~j) . 

VI. MSUGRA 

Before we proceed to present our results, let us first mention the values of the parameters 
used in our calculation of B — > <pK decay amplitudes in SUGRA models since the BRs 
depend sensitively on them. We use pa,h = 1, <Pa,h = —68°, m s (2GeV) = 122.5 MeV and 
a CKM fit giving sin 2(3 = 0.73 and 7 = 59°. (If we increase 7, the BR decreases, e.g. 
for 7 = 79°, the SM BR decreases by ~ 2%.) The SM BR is 4.72 x 10~ 6 and the weak 
annihilation contribution is small (~ 10 %). 

In Table HJ we give the numerical results for two different values of the mSUGRA pa- 
rameter tan/3 cases i.e. tan/5 = 40 and 10 and for different m^j and A Q in the small 
weak annihilation case. For simplicity, we set oa = ^ in the calculation of Table I. We use 
large phases for other parameters but still satisfy the EDM constrains. For example, for 
m 1/2 = 400 GeV and A = 800 GeV with tan/3 = 40, we find that fa = 70°, (f) 3 = 33° and 
M = —13° (at the weak scale) satisfy the EDM constraints (for reasons discussed in detail 
in j^l). The phase a a has a very small effect on S^Ks an d this effect becomes smaller as the 
magnitude of A decreases. Thus a different value of a a, in the above fit, can change S^Ks 
by ±4% for Aq = 800 GeV. This change is even smaller for smaller Aq, e.g. for Aq = 200 
GeV, the change in S^Ks is l ess than 2%. So far we have not specified m for our results. 
The values of mo, for different and A , are chosen such that the relic density constraint 
is satisfied. We also satisfy the Br [b — > s + 7] constraint and the Higgs mass constraint. The 
S(f>K s values shown in the table are the minimum that can be reached for given mi/2 and Aq. 

It can be seen from Table I that the S^Ks values in mSUGRA differ only slightly from the 
SM prediction which is sin 2/3 evaluated using just the CKM phase. The branching ratios of 
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B — > (f)K decays also do not differ much from the SM prediction. Even if one went to the 
large weak annihilation region to accommodate the large branching ratios, S^Ks would still 
be similar to the numbers in Table I. Therefore, mSUGRA can not explain the large BR 
and the 2.7a difference between the S^k s an d the Sj/^k s experimental results. The reason 
is that, in mSUGRA, the only flavor violating source is in the CKM matrix, which cannot 
provide enough flavor violation needed for the b — >• s transition in B — > <pK decays. In the 
next section, we will search for the minimal extension of mSUGRA that can solve both the 
BR and CP problems of B — > <pK decays. 



tan/? 


10 


40 


\M 


800 


600 


400 





800 


600 


400 





m i/2 = 400 


0.74 


0.74 


0.73 


0.73 


0.71 


0.70 


0.70 


0.69 


m\j2 = 500 


0.74 


0.74 


0.74 


0.74 


0.72 


0.72 


0.72 


0.71 



TABLE I: S^Ks with small weak annihilation at tan/3 = 10 and 40 in mSUGRA. The values are 
the minimum that can be reached subject to all other experimental constraints 



VII. SUGRA MODEL WITH NON-UNIVERSAL A TERMS 

In the last section, we showed that mSUGRA contributions to B — > <j)K decays are 
negligible and thus mSUGRA needs to be extended if it is to explain the experimental 
results of B — > <pK decays. It is obvious that some non-universal soft breaking terms which 
can contribute to the b — > s transition are necessary. There are two ways of enhancing the 
mixing between the second and the third generation: one can have non universal terms in 
the squark mass matrices or in the A U,D matrices of Eq. (jl4j) . However, in a GUT model, at 
least the Standard Model gauge group must hold at Mq and hence the only squark m| 3 that 
can occur is either left-left or right-right coupling. As discussed in j^j], such non-universal 
terms produce only small effects onB-> <pK decays. Thus we are led to models with left- 
right mixing which can occur in the A U,D matrices as the simplest possible non-universal 
term relevant to B — > <pK decays. In this work then, we choose a model with non-zero (2,3) 
elements in the trilinear coupling A terms of Eq. ()14j) to enhance the left-right mixing of the 
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second and the third generation. The A terms with non-zero 23 elements can be written as 

A U,D = A Y U ' D + AA U ' D (23) 

where AA U > D are 3 x 3 complex matrices and AA^f = \AAfj D \e^ D . When AA U ' D = 0, 
mSUGRA is recovered. For simplicity, we will discuss first the case of non-zero AAg and 
non-zero AAg for tan/3 = 10 and 40. In both cases, all other entries in AA U,D are set 
to zero. The other parameters are same as in the mSUGRA case. We will set the phases 
such that the EDM constraints are obeyed. At the GUT scale, we use a diagonal Yukawa 
texture for Y u , while Y D is constructed as VY® where V is the CKM matrix and Y® 
is the diagonalized matrix of the down type Yukawa. The phenomenological requirements 
for the Yukawa matrices are that they produce the correct quark masses and the correct 
CKM matrix. Any other Yukawa texture which satisfies the same requirements can be 
obtained through unitary rotations. Therefore, our results can be recovered with other 
Yukawa textures if our A terms are rotated along with the Yukawas. 

In the calculations of decay amplitudes, we will use QCD parameters for the small weak 
annihilation re gion (see the last section) where the theory is reliable. In general it is pos- 
sible that (see |2jj for the calculational details of weak annihilation) the new physics can 
change the behavior of annihilation contributions when the relevant Wilson coefficients can 
be reduced or increased significantly. However, in our case with non-zero AA\f Z2 terms, the 
SUSY contribution mainly affects the Wilson coefficients Cg g (7 7 ) (possibly also C8 S (7 7 )) and 
these coefficients will not change the annihilation contributions compared to what we have 
in the SM calculation and thus our previous conclusion about the annihilation terms still 
holds. 

Case I: \AAg\ = \AAg\ and <f>g ^ <f>g 

We show our results for \AAg\ = \AAg\ but 4>g ^ <j)g with tan/? = 10 in Table HU 
We note that |Ay4^ ( - 32 ^| is an increasing function of m^j. The phases (ftg and <f)g are 
approximately —30° and (75 ~ 115)°, respectively. The other SUSY phases are: 0i ~ 22°, 
</>3 ~ 31° and <f)^ ~ —11°. In addition, as mentioned above, the phase of A , i.e. a^, is set 
to be 7T. 

The Br[.B~ — > (j)K~\ is ~ 10 x 10~ 6 in the parameter space of Table II. We satisfy 
all other experiment constraints. We see that SUGRA models can explain the large BR 
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\A \ 


800 


600 


400 





\AA D 1 


mi/2 = 300 


-0.50 


-0.49 


-0.47 


-0.43 


~ 50 


m i/2 = 400 


-0.43 


-0.40 


-0.38 


-0.36 


~ 110 


m i/2 = 500 


-0.46 


-0.46 


-0.44 


-0.31 


~ 200 


m i/2 = 600 


-0.15 


-0.13 


-0.04 


0.05 


~ 280 



TABLE II: S$ Ks at tan/? = 10 with non-zero ^23 an d ^32 
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600 


400 





m i/2 = 300 


1.2% 


-3.7% 


1.4% 


-3.6% 


1.7% 


-3.6% 


2.2% 


-3.5% 


m i/2 = 400 


1.9% 


-3.5% 


2.0% 


-3.4% 


2.2% 


-3.3% 


2.3% 


-3.3% 


m i/2 = 500 


2.6% 


-3.5% 


2.6% 


-3.6% 


2.5% 


-3.5% 


2.4% 


-3.2% 


m i/2 = 600 


2.0% 


-2.8% 


2.1% 


-2.7% 


2.1% 


-2.5% 


2.2% 


-2.2% 



TABLE III: A-» s + 7 (left) and A$ K ^ (right) at tan/? = 10 with non-zero ^23 and A^- 

and S(f,K s of the B — > <j)K decay modes even in the small annihilation region. Comparing 
with Eq.(JIJ, one sees that the values of S^Ks i n the Table are within la range of the 
experimental measurement. Reducing the Br[B~ — > (j)K~\ allows one to lower S ( j > K s - For 
example, for A = and mu2 = 600 GeV, by adjusting (f)® 2 i S(/>k s can be reduced to -0.05 
with Br[B~ — > <pK~] ~ 9 x 10 -6 . In Table ITTT1 we show the direct CP asymmetries of the 
B~ — > (f)K~ decay, i.e. A^k^i using the same parameters as Table ITTl The CP asymmetry is 
around — (2 ~ 3)% and agrees with the experimental result shown in Eq.©. This prediction 
depends on the choice of 4>a,h in Eq. (j2*T|) . For example, if we use 4>a,h = 28°, we generate a 
large A^k^ ~ 27%. We find that there exists a reasonably large range of 4>a,h where we can 
satisfy the current bound on A^k^- For example, at m^j = 300 GeV and A = 800 GeV 
where the SUSY contribution is the largest, we find that A^k^ varies from 9% to —4% when 
<t>A,H varies from —100° to —50° (for simplicity, we set 4>a = 4>h)- In addition, since the 
annihilation contribution is small in that range, we find that the branching ratio is around 
(9.5 ~ 11) x 10~ 6 . The CP asymmetry of b — > 37 is ~l-3%. The present experimental errors 
for C^Ks are still large. For this model, C$k s ~ — A^kt, which may be tested by future 
data. 
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-0.40 
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-0.38 


10.0 


-0.33 
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10.0 


mi/2 = 400 


-0.11 


8.0 


-0.05 


8.0 


0.04 


7.9 


0.28 


8.0 


mi/2 = 500 


0.07 


6.0 


0.16 


6.1 


0.24 


6.1 


0.37 


6.2 


m i/2 = 600 


0.37 


6.2 


0.44 


6.2 


0.49 


6.2 


0.58 


6.2 



TABLE IV: S (j)Ks (left) and Bt[B~ -> <j>K-\ x 10 6 (right) at tan j3 = 40 with non-zero AA% 3 and 
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-5.6% 


-3.4% 


-5.2% 


-3.3% 


-3.6% 


-2.6% 


mi/2 = 400 


-3.0% 


-3.0% 


-2.1% 


-2.9% 


-1.7% 


-2.6% 


-0.7% 


-1.8% 


m i/2 = 500 


-0.5% 


-2.9% 


-0.4% 


-2.5% 


-0.2% 


-2.2% 


0.2% 


-1.7% 


mi/2 = 600 


0.2% 


-1.7% 


0.3% 


-1.4% 


0.4% 


-1.2% 


0.6% 


-0.8% 



TABLE V: A h -^ s +y (left) and A^k^ (right) at tan f3 — 40 with non-zero A 23 and A 



32- 



In Table IIVI and Table we give our results for tan (5 — 40 with non-zero AA^^) ■ 
The phases <j>g and <f>& are -(70 ~ 0)° and (80 ~ 110)°, respectively, fa ~ (25 ~ 60)°, 
03 ~ 25° and <^ ~ -8°. The off-diagonal elements \ AAg {32) \ vary from 90 GeV to 250 GeV 
as mi/2 increases. We compare Table ITVl with the results for tan/3 = 10 shown in Table 
|n]and we see that only low mi/2 can satisfy experimental data for tan/3 = 40. The most 
important reason for this is that left-right mixing of the second and the third generation 
decreases significantly with increasing tan (3. This comes about as follows. The RGE running 
of ^213(32) i s n °t sensitive to tan /?. Therefore, for the same size of ^23(32) m P u t at the GUT 
scale, the weak scale values of do not differ much for different tan (3. However, the 

A D term enters into the down squark matrix after electroweak symmetry breaking when 
ifi (see Eq.(|14j)) grows a vacuum expectation value proportional to cos/3. Hence the left- 
right mixing between the second and the third generation in the down squark matrix will 
be smaller for large tan/3. For low my2 this reduction can be compensated by increasing 
the magnitude of Ag^y For example, for rai/2 = 300 GeV, we use | ^.23(32) I ~ 90 GeV 
in this case compared to 50 GeV in the case of tan/3 = 10 (see Table ITT|) . The chargino 
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diagram contribution increases with tan/5 and can help to generate a large BR. But for 
large m^, when the chargino contribution goes down, 1^4.23(32)1 niust become much larger. 
However, as A 2 ^2) increases, the pseudoscalar Higgs mass becomes small at the same time 
(but fj, does not get smaller), which prevents I ^23(32) I f rom having an unlimited increase. For 
example, for mi/ 2 = 600 GeV and A = 800 GeV, I ^-23(32) I = 250 GeV generates wla = 580 
GeV which is still allowed for the dark matter constraint to be satisfied in the f <-> x° co- 
annihilation channel. If I ^-23(32) I * s increased more, the pseudoscalar mass gets smaller and 
the dark matter constraint can still be satisfied due to the available X1X1 ~* A —* f f channel. 
But with a further reduction of the pseudoscalar mass by increasing | ^23(32) I further, this 
channel goes away when rriA < 2m^o and we must again satisfy the relic density using the 
stau-neutralino co-annihilation channel. However, the improvement of S<f,K s in this scenario 
is small. For example, for the point mentioned above, I ^23(32) I can ^ e increased to around 
480 GeV with relic density in the f <-> x° channel but S^k s can only be reduced from 0.37 
(see Table HVjl to 0.22 with the same branching ratio. Thus, the S^,k 8 and the branching 
ratio still cannot be satisfied. 

If we use (P23 — 032 (equal phases) we have one less parameter, but that choice will 
not be able to satisfy experimental results. The reason is that the weak phase from the 
gluino contributions in the Wilson coefficients C 8g and the weak phase from C$ g will cancel 
when <f>2 3 = (j)§2 because C 8g depends on A 23 but C 8g depends on (Af 2 )*. For example, for 
tan/? = 10 we find that S^Ks ~ 0-7 since the gluino contribution dominates at lower tan/?. 
At tan j3 = 40, S 1 ^ can reach 0.45 since the chargino contribution is larger at higher tan j3, 
but this is not enough to satisfy the data. 

Case II: \AAg\ = \AA&\ and cg z = <pg 

In this section we discuss the case AA 23 ^ = but AA 23 ^ 7^ 0. The phases used are 
similar to those used in first two cases except 023 = 032- This case is more complicated than 
the ^3(32) 7^ case. We find that it is easier to start by comparing them. 

The first important change is that the AA 32 contribution is much smaller than the AA 32 
contribution to the mixing between the second and the third generation in the down squark 
mass matrix due to the suppression by the second generation Yukawa coupling in the RGE 
of A 32 . (Thus our choice of 0^ = <p 32 has no loss of generality.) Consequently, the size 
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600 


400 





1^(32)1 (GeV) 


mi/2 = 300 


0.03 


8.4 


0.04 


9.0 


0.01 


8.0 


0.17 


8.0 


~ 300 


m x /2 = 400 


-0.07 


8.5 


-0.03 


8.4 





7.1 


0.32 


6.3 


~ 600 


m \/2 = 500 





6.5 


0.07 


6.4 


0.18 


6.0 


0.44 


6.1 


~ 800 


m i/2 = 600 


0.27 


6.1 


0.30 


6.1 


0.35 


6.1 


0.51 


5.9 


~ 1000 



TABLE VI: S (j)Ks (left) and Br[B~ -> ^>K~] x 10 6 (right) at tan/3 = 40 with non-zero AA% 3 and 





800 


600 


400 





1^(32)1 ( GeV ) 


m i/2 = 300 


0.17 


6.5 


0.16 


6.3 


0.32 


6.1 


0.60 


5.2 


~ 300 


m 1 / 2 = 400 


0.37 


4.7 


0.39 


4.6 


0.46 


4.3 


0.62 


4.3 


~ 550 



TABLE VII: Sqks (left) and Bt[B~ -» <j)K~] x 10 6 (right) at tan/3 = 10 with non-zero AA% 3 and 



of the Wilson coefficient C$ g is significantly reduced. Although AA^ 3 still contributes, that 
contribution is also reduced (compared to AA^) due to the RGE. Therefore, compared with 
the first case the total SUSY contributions are reduced especially for tan (3 = 10 and thus it 
becomes harder to fit the experimental results. 

Another important change is the roles of some experimental constraints which are not 
important in the first case in the sense that they do not prevent the SUSY contributions 
from increasing, or at least their limits are not reached when we have solutions satisfying 
the B-decay data. Below are some comments concerning this: 

1. For tan/3 = 40 and low m^, i.e. 300 GeV, the Br [5 — > Xg'j] will constrain the 

l 23(32) 



size of AAwfmy, . This is why the S^Ks an d the branching ratio fit is not as good as the 



corresponding one shown in Table |H] for the A!^^) 7^ case. 

2. When m 1( / 2 increases, the size of AA^^) also needs to be increased. But three other 
additional constraints are present, i.e AM K and ex from the K° — K° mixing and the mass 



of smallest up squarks (right-handed stop) nif . For example, for rr 
(and rriQ = 431 GeV by the relic density constraint) we get rrig ~ m f - 



1/2 



is the average squark mass and m~ g is the mass of the gluino, see 
we find that 



500 and A = 600 
1000 GeV (where m- q 



39j for more details) and 



Re(5f 2 )! L | = 7.1 x 10" 2 and J\lm{6f 2 ) 2 LL \ = 9.7 x 10~ 3 which are allowed by 
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the experimental bounds on AMk and ex |39j (the sizes of {Sf 2 )LR, ($12) RL an d (^12) RR are 
around 10~ 8 ~ 10~ 7 and thus these bounds can be safely ignored in our case). 
3. The situation for the right-handed stop mass is similar to the pseudoscalar Higgs case 
we mentioned at the end of Case I. We use the f <-> x° channel to satisfy the dark matter 
constraints. Although it's possible to use a larger A 2 ^ 2 ) which consequently reduces 
more and then opens the tn <-» x° channel, the room is small due to the smallness of m^o. In 
addition, the Mk and the ex bounds become harder to satisfy when rrtf is small. Therefore, 
as in the case where the pseudoscalar Higgs mass becomes small, possible improvements can 
not satisfy the experimental results of both S<j,K s an d Br[B~ — > <pK~}. 

A further difference is the tan/3 dependence. In Case I, as was discussed above, gluino 
contributions depend inversely on tan/5 due to the way that A^^) en t ers hito the down 
squark mass matrix. But in this case, the gluino contributions are reduced significantly and 
the chargino plays a more important role, which will be enhanced by tan/3. Therefore, in 
this case, we see that larger tan/3 can satisfy the experimental results at small mm, but 
small tan (3 cannot and that is why we have only given results in Table IV11I for two values 
of mi/2 at tan/3 = 10 since higher m-1/2 cannot improve the situation. 

We also comment concerning Aq and its phase. So far, we have used the phase ir for Aq. 
We find that using a different phase will not improve the results greatly. In general, the 
improvement is at a few percent level. (This holds also for case I.) For example, in Case II, 
for tan/3 = 40, mi/ 2 = 400 and A = 800, we find that using a a ~ —95° can improve S$k s 
from -0.04 to -0.06. 

Finally we note that the values of Ab^ s + 1 and A^k* remain small, i.e. Ab^ s+1 and A^k* 
are -(3 ~ 0)% and -(3 ~ 1)% at tan/3 = 10, and -(5 ~ 0)% and -(3 ~ 1)% at tan/3 = 40. 

VIII. CONCLUSION 

In this paper, we have probed the B — > (j)K decays in SUGRA models with CP violating 
phases to explain the discrepancy between the experimental measurements and the SM 
predictions of the CP asymmetry of B° — > <pKs and the branching ratios of the B — > (j)K 
decays. We have calculated the CP asymmetries of B~ — > <pK~ and B — > X s ^. In our 
analysis, we implemented all relevant experimental constraints, e.g. Br [B — > X s j], relic 
density, K° — K° mixing parameters and electron and neutron EDMs. We used the improved 
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QCD factorization method |26|, |2jj for the calculation of decay amplitudes. 

The SM not only can not explain the CP asymmetry of B° — ► 4>Ks, it also fails to satisfy 
the Br [B — > <pK] data barring the region of large weak annihilation where the theory is most 
suspect. We then studied the mSUGRA model and found that it also has the same problem. 
Therefore, if the current experimental results continue to hold in the future, it will signal the 
first significant breakdown of the Standard Model and also of mSUGRA. This conclusion 
is important in the sense that one needs to construct a more complicated SUGRA model 
to satisfy experimental measurements which will provide important guidance to our future 
research on SUSY models and their signals at the accelerator experiments. 

We have considered the extension of the mSUGRA model by adding non-universal A 
terms. For a GUT theory, the only natural choice is to have a left-right mixing between 
the second and the third generation in the up or down quark sectors i.e. AA^ 13 and AA^ 13 
terms. We have examined thoroughly several different possibilities in this extension and 
their theoretical predictions and have found a large region of parameter space where all 
experimental results can be satisfied, including the CP asymmetries and branching ratios of 
the B — > cf)K decays. This result is obtained without resorting to large weak annihilation 
amplitudes and so is based on reliable calculations of hadronic decays, and thus provides 
useful hints for the study of hadronic B decays. Further the size of AA^ needed is the 
same as the other soft breaking terms, and so is not anomalously small or large. Thus, this 
study also can provide important phenomenological information not only for accelerator 
physics but also for building models at the GUT scale and for exploring physics beyond 
it. In this connection, models based on Horava-Witten M-theory can naturally give rise 



to non-zero values of AA 23 . In previous work |42| it was shown that it was possible to 
construct a three generation model with SU(5) symmetry using a non-standard embedding 
based on a torus fibered Calabi-Yau three fold with a del Pezzo base dPj. The model 
allowed Wilson line breaking to the Standard Model at Mq, and also had vanishing instanton 
charges on the physical orbifold plane. If in addition one assumed that the 5-branes in the 
bulk clustered around the distant plane, one could explain without undue fine tuning the 
general structure of the quark and lepton mass hierarchies and obtain the LMA solution for 
neutrino oscillations j^, 143I [j^ . One can show that the model naturally gives rise to AA 23 
non-universal terms as required by the current B-factory data, while keeping the squark 



mass matrices essentially universal. This model will be discussed elsewhere [45]. 
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